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Abstract 

The direct detection of dark matter through its elastic scattering off nucleons is among the most 
promising methods for establishing the particle identity of dark matter. The current bound on the 
spin-independent scattering cross section is <r SI < 10 zb for dark matter masses m x ~ 100 GeV, with 
improved sensitivities expected soon. We examine the implications of this progress for neutralino 
dark matter. We work in a supersymmetric framework well-suited to dark matter studies that is 
simple and transparent, with models defined in terms of four weak-scale parameters. We first show 
that robust constraints on electric dipole moments motivate large sfermion masses rh > 1 TeV, 
effectively decoupling squarks and sleptons from neutralino dark matter phenomenology. In this 
case, we find characteristic cross sections in the narrow range 1 zb < ct si < 40 zb for m x > 70 GeV. 
As sfermion masses are lowered to near their experimental limit rh ~ 400 GeV, the upper and 
lower limits of this range are extended, but only by factors of around two, and the lower limit 
is not significantly altered by relaxing many particle physics assumptions, varying the strange 
quark content of the nucleon, including the effects of galactic small-scale structure, or assuming 
other components of dark matter. Experiments are therefore rapidly entering the heart of dark 
matter-favored supersymmetry parameter space. If no signal is seen, supersymmetric models must 
contain some level of fine-tuning, and we identify and analyze several possibilities. Barring large 
cancellations, however, in a large and generic class of models, if thermal relic neutralinos are a 
significant component of dark matter, experiments will discover them as they probe down to the 
zeptobarn scale. 

PACS numbers: 95.35.+d, 12.60.Jv, 14.80.Nb 
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I. INTRODUCTION 



Now is an exciting time for the study of non-baryonic cold dark matter (DM). The exis- 
tence of DM offers a consistent explanation for a wide variety of astrophysical observations 
on disparate length scales. At the same time, there are many viable dark matter candi- 
dates, motivating a diverse and active research program searching for the direct detection 
of dark matter scattering off nuclei, the indirect detection of dark matter annihilation or 
decay products, dark matter production at colliders, and the effects of dark matter particle 
properties on structure formation and other astrophysical phenomena [l|-[3|. 

In this work, we focus on direct detection. We will discuss and present results for both 
spin-independent and spin-dependent scattering, but will focus particularly on probes of 
the spin- independent dark matter- nucleon scattering cross section a sl . Although no defini- 
tive signal of dark matter has been established so far, recent results from the DAMA jiJ], 
XENON10 Q, CDMS-II @, CoGeNT [7], and XENON100 |M3 experiments have sparked 
great interest. Moreover, significant progress for spin- independent cross sections is expected 
from several sources, including XENON100 (Tj, CRESST Q, LUX [13], XMASS Q, 
mini-CLEAN/DEAP [l5|, WArP 0, and superCDMS |T3|. Current limits are so stringent 
that there is no agreed upon unit convention for expressing them. We advocate the use of 
zeptobarns (1 zb = 10~ 9 pb = 10~ 45 cm 2 ), which is more compact and no less transparent 
than the alternatives. The current limit is cr SI < 10 zb for 100 GeV dark matter 



In the near future, results from 100 kg-scale detectors may improve this sensitivity by an 
order of magnitude to the zeptobarn scale. Future ton-scale detectors are being planned 
to extend sensitivities even beyond that, with the ultimate goal of probing yoctobarn cross 
sections, where experiments will reach the background-free limit for non-directional detec- 



tors imposed by irreducible neutrino backgrounds [la . 1 1 91 ] (and also the outer limits of the 
metric prefix system). 

The single most studied dark matter candidate is the lightest neutralino x predicted by 
supersymmetry (SUSY) [20-24]. SUSY provides a solution to the gauge hierarchy problem 
and gauge coupling unification, and with .R-parity conservation, the lightest supersymmetric 
particle (LSP) is stable. The LSP in SUSY theories is often the lightest neutralino, a mixture 
of the superpartners of the neutral electroweak gauge and Higgs bosons. It is an excellent 
DM candidate, not only because SUSY is a well-motivated possibility for new physics, but 
also because the lightest neutralino is a weakly-interacting massive particle (WIMP), and 
so naturally has the appropriate thermal relic density. 

What are the possible values of <r SI for neutralino dark matter? In full generality, the al- 
lowed range varies from very large values that are far above current experimental bounds 25- 



271 ] to very small values that are far below the sensitivity of foreseeable experiments [28]. At 



the same time, there is some sense in which the extreme values occur only in special regions 
of SUSY parameter space. In this study, we aim to answer a more subtle (and less well- 
defined) question: what are the characteristic values of a for neutralinos? As experiments 
improve, is every factor-of-ten improvement in cx SI equally important, or are there certain 
thresholds that should be considered particularly significant? In other words, what is a in 
the heart of SUSY parameter space? 

To answer these questions, one must define preferred models, reducing the general param- 
eter space to some subspace by appealing to some organizational principle. One approach 
is to work with a low-energy effective theory that includes the dark matter particle but not 
other supersymmetric degrees of freedom. Dark matter-quark couplings are included with, 
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for example, four-fermion interactions suppressed by some large mass scale. This 

approach has the advantage of being model-independent, with results applicable to theo- 
ries beyond SUSY, and has been adopted in a number of recent studies motivated by the 



possibility of light (GeV-scale) WIMPs [29H32| . The effective theory approach is, however, 



inapplicable to the generic SUSY theories of interest to us, in which the dark matter mass 
is not hierarchically smaller than those of the other superparticles. Also, effective theory 
studies typically assume that the same operators govern both annihilation and scattering, 
and so neglect the annihilation to leptons and gauge bosons that may be important in SUSY 
theories. 

Another approach has been to focus on a particular high-energy framework, incorporating 
constraints from a battery of known observables, and scanning over the reduced, but still high 
dimensional, parameter s pac e. The most popular framework for these analyses is minimal 



supergravity (mSUGRA) [33M37l| , with its cosmologicalfy-preferred regions, including the co- 



annihilation region at low universal scalar mass mn 38l. l39l| . the focus point region at high 



mo |40H42| . and the A-funnel region at high tan/3 ]39|,j43j. Much work has been devoted to 
this framework; for recent analyses, see, e.g., Refs. [44j-|47|. This approach has the advantage 
that there is a manifestly consistent high-energy formulation throughout parameter space, 
but there are also drawbacks. The relation of the fundamental high-energy parameters 
to physical observables, such as superpartner masses, is typically opaque, at least to non- 
experts. In addition, the framework allows one to relate dark matter properties to b — > sj, 
the anomalous magnetic moment of the muon, and other a priori unrelated observables, and 
it is not clear how robust these connections should be. Finally, not all of the assumptions of 
mSUGRA are well-motivated, and it is not always easy to determine how dark matter results 
would change if some of the assumptions were relaxed or another framework adopted. For 
example, are the promising predictions for a sl in focus point supersymmetry necessarily tied 
to large scalar masses, with their negative implications for colliders? Are there consequences 
of the connection between sfermion and Higgs masses in mSUGRA that fail to translate to 
more general theories? 

In this work we take a complementary approach. Although we will apply some constraints 
and characteristics from high-energy motivations, we work in a manifestly low-energy frame- 
work, with models defined by a few weak-scale SUSY parameters that are easily related to 
physical observables, typically the physical masses of certain superpartners. The model 
framework is particularly well-suited to dark matter studies for the reasons given below, 
and although specific to SUSY, our results will apply to a broad class of SUSY models. In 
addition, it will allow us to explore the entire parameter space and understand the physics 
behind the results in a detailed and relatively transparent way. This in some ways resembles 



the approach taken in the so-called phenomenological MSSM (pMSSM) [48145 though our 
exclusive interest in neutralino dark matter allows for a more directed study of a smaller pa- 
rameter space. We will find that many correlations found in mSUGRA analyses are broken 
in these more phenomenological models. In particular, the large signals for direct detec- 
tion found in focus point supersymmetry are generic and not necessarily tied to large scalar 
masses, with highly favorable implications for the immediate future of neutralino direct 
detection. 

This paper is organized as follows. In Sec. [TT1 we define our model framework and discuss 
its input parameters and assumptions. We also discuss the most relevant constraints, notably 
those from bounds on the electric dipole moments (EDMs) of the electron and neutron. In 
Sec. II I II we analyze the dependence of the thermal relic density on the model parameters, 
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and in Sec. IIVI we briefly review the form of cx SI and explain its expected behavior in our 
model framework. In Sec. |V]we present our main results for the characteristic values of <r SI , 
the possibility of fine-tuned cancellations, and the implications for dark matter discovery. 
In Sec. ED we explore the effects of relaxing various particle physics, nuclear physics, and 
astrophysical assumptions, and we discuss the prospects for spin- dependent scattering in 
Sec. IVHI We summarize our conclusions in Sec. IVIIII 



II. MODEL FRAMEWORK 

A. Input Parameters and Assumptions 

We consider SUSY models with minimal field content that are specified by five parameters 
defined at the weak scale ~ 1 TeV: 

rh, Mi, (jl, rriA, tan/3 . (1) 

In order, these are the unified sfermion soft SUSY-breaking mass, the Bino mass, the Hig- 
gsino mass, the pseudoscalar Higgs boson mass, and the ratio of Higgs boson expectation 
values. 

The rest of the model is specified by the following three assumptions: 

1. Gaugino mass unification. We assume that the gaugino masses are unified at the grand 
unified theory (GUT) scale, implying the weak-scale gaugino mass relation 

Mi : M 2 : M 3 ~ 1 : 2 : 6 . (2) 

2. No left-right mixing. We assume there is no significant left-right sfermion mixing, 
which fixes the "A-term" tri-linear scalar couplings. For example, the top squark 
mass matrix simplifies as follows: 



m| + mf m t (A t - /i cot 0) \ / m 2 + m 2 
m t (A t — /j, cot f3) ml + ntf J \ 



2 

(3) 




3. Correct electroweak symmetry breaking. This fixes the soft Higgs scalar mass param- 
eters m 2 H d and m\ 2 through the (tree-level) conditions 

2 _ m H d ~ m H u tan2 P + (1 ~ tan 2 /3) 
^ ~ tan 2 /3-l W 

m l = ™ 2 H d + m 2 H u + 2/i 2 (5) 

m 12 = m A sm P cos P ■ (6) 

The model framework defined above is simple, but well-suited to studies of neutralino 
dark matter. With the assumption of gaugino mass unification, the parameters Mi, \x, and 
tan (3 completely determine the neutralino mass matrix, and so are natural input parameters 
for this study. These also fix the chargino mass matrix and the gluino mass. The remaining 
input parameters m and rriA determine the masses of all the remaining particles. In our 
numerical work, we use SOFTSUSY 3.1.2 [52| to generate the superpartner spectrum and 



micrOMEGAS 2.2 [53|, |54[ for relic density and direct detection calculations. 
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We conclude this subsection with a few further comments about the model parameter 
space and assumptions. Throughout this analysis, we will assume that the neutralino freezes 
out with the correct thermal relic density Q x = Qbm — 0.23. Our results will be remarkably 
insensitive to this assumption, as discussed in Sec. IVIF} but this constraint further simplifies 
the parameter space. As discussed in Sec. 1111} for a fixed m x , the relic density determines /i 
up to a discrete choice, which may be taken to be the sign of /i. The input parameters then 
become 

fh, Mi, \i, rriA, tan/3 nx ~°' 23 ) ^ sign(/i), m^, tan/3 . (7) 



We will typically plot quantities as functions of m x for the two choices of sign(/i) and fixed 
rh, rriA, and tan/3. 

The assumption of gaugino mass unification is motivated by the unification of gauge 
couplings and forces, a key virtue of supersymmetry. It implies that the LSP does not have 
a large Wino component, which restricts the generality of this analysis. On the other hand, 
a neutralino with a significant Wino component typically annihilates too efficiently to have 
the correct thermal relic density, and so our assumption that the neutralino is a thermal 
relic already disfavors this possibility [55| . Moreover, while neutralino dark matter with non- 



unified gaugino masses has been discussed [56|, |57( , and more complicated GUT-breaking 



schemes are possible jlsl 591, the gaugino mass unification scenario used here remains a 
standard assumption for many high-scale models. However, it is worth noting that our 
results do not require strict gaugino mass unification, but only that M2 and M3 are at least 
somewhat larger than Mi. 

The assumptions of a unified weak-scale sfermion mass fh and negligible left-right mixing 
are not motivated by high-energy theories of SUSY-breaking; instead, they are motivated 
by the relative insensitivity of dark matter phenomenology to these assumptions. We will 
explore the implications of deviations from the above assumptions in Sees. IVI Bl and IVI CI 
respectively. We note here, however, that the dependence on left-right mixing is small, 
and the neutralino interactions of greatest interest here depend primarily on the squarks 
in most regions of parameter space q In many theories the assumption of a unified squark 
mass is reasonably accurate, and so this parametrization characterizes these models well, 
with fh playing the role of the characteristic squark mass. Throughout this study, we 
neglect sfermion flavor mixing and one-loop corrections to sfermion masses, as these are 
sub-dominant effects in neutralino interactions. 

Finally, in our model framework, correct electroweak symmetry breaking is achieved at the 
expense of unifying the soft Higgs scalar masses with those of the squarks and sleptons. The 
unification of the Higgs scalar masses with the soft squark and slepton masses is assumed in 
mSUGRA, but it is one of the least motivated aspects of that framework. Flavor constraints 
motivate flavor-independent squark and slepton masses, but do not constrain the Higgs 
scalar masses. Even in GUTs, the Higgs scalars and sfermions are typically in different 
multiplets. Our assumption regarding the Higgs mass parameters follows other analyses of 
non- unified Higgs models (NUHM) [60 631 . 



The relic density only depends strongly on slepton masses in the presence of slepton co- annihilation or 
very light slepton masses. 
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B. Constraints 



1. Flavor and CP Violation 



Theories that introduce new particles at the weak scale run the risk of violating stringent 
bounds on low-energy flavor- and CP-violating observables. For SUSY, there are mediation 
mechanisms that lead to generation-independent slepton and squark masses, which allow 
these theories to satisfy constraints on flavor violation. These will also be satisfied in our 
models, given the assumption of a unified weak-scale sfermion mass fh. 

Even in these cases, however, constraints from EDMs, which conserve flavor but violate 
CP, are not necessarily satisfied. In a general SUSY model, the EDMs of the electron and 
neutron receive contributions from several sources. Left-right sfermion mixing can lead to CP 
violation, but this is absent in our model framework. However, there are also contributions 
from penguin diagrams with gauginos, Higgsinos, and a sfermion of fixed chirality in the 
loop, which are proportional to sin^cp, where 0cp is the mismatch in phases of the gaugino 
masses and \i parameter. In the gauge eigenstate basis, the dominant diagrams are those 
with neutral Winos, neutral Higgsinos, and left-handed sfermions in the loop, and a similar 
one with charged Winos, charged Higgsinos, and sneutrinos in the loop. The first diagram 



leads to the EDM contribution 64 



-era/j 2 2 \M1\1\ tan/3 sin </>cp K^{m^ \fi 



\ \Mo\ 2 ) 



(8) 



7i -"m- 

JL 



where Kn is a kinematic function that may be approximated as Kn ~ 1/(16 

The current bounds on the electron and neutron EDMs are d e < 1.6 x 10~ 27 e cm [6f 
and d n < 2.9 x 10~ 26 e cm j67[. Assuming m u = 3 MeV, = 5 MeV, the naive quark 
model relation d n = (Ad^ — d u )/3, and neglecting cancellations between different diagrams, 
we find 



\M 2 fi\ 

m k 
\MzfA 

2 

m- 



6.5 TeV 



2.1 TeV 



m 



5.5 TeV 



m d L 



tan/3 
10 

tan/3 
10 

tan j3 
10 



sin (p C p 



sin (p CP 



sin (p CP 



<1 



<1 



<1 . 



(9) 
(10) 

(11) 



In our model context, with unified scalar masses and M 2 ~ 2Mi, these imply 

n 2 



\Mxll\ 

rh 2 



9.2 TeV 



m 



tan /3 

To" 



sin (pep 



<1 . 



(12) 



Note, however, that even without assuming unified slepton and squark masses, the con- 
straints on mg L and are similar, and the constraint on is only somewhat less 
restrictive. 

This is a robust and numerically stringent constraint. Even for the low values \Mi\, ~ 
100 GeV and moderate tan/3 ~ 10, this implies [1 TeV/m] sin0cp ~ 1- Without a mecha- 
nism for aligning the CP- violating phases of different SUSY parameters, we naturally expect 
shi(/>cp ~ 1, implying m > 1 TeV. Lower values of m are possible, of course, but require 



6 



small values of 0cp- As we will see, for m ^ 1 TeV, the sfermions are largely decoupled 
from neutralino interactions, leading to many simplifications in the analysis of neutralino 
annihilation and scattering. We will present results for a variety of m up to 2 TeV. For this 
largest value, arbitrary CP-violating phases are allowed for \M\\, |/i| ~ 400 GeV, and only 
moderate suppression from sin0 CP < 0.2 is necessary for |Mi|, < 1 TeV. 

There are several other constraints on flavor and CP violation that are often used to eval- 
uate the viability of SUSY theories. The most notable are the constraints on the anomalous 
magnetic moment of the muon a M = (g M — 2)/2 [6$] and the rare decays b — > s^y [691 and 



b — > [70]. We will neglect these constraints in our analysis, as large sfermion masses 

motivated by EDMs largely avoid these constraints. The experimental value for B(b — > 37) 
is consistent with the SM prediction (71, 72], and B(b — >■ sfi + [i~) has only an upper limit. In 
both cases, large sfermion masses suppress the SUSY contributions to these observables suf- 
ficiently to avoid the constraints. The contribution to is also suppressed for large sfermion 
masses, implying that models in our framework do not alleviate the a M discrepancy [73] , but 
do not exacerbate it either. 



2. Higgs Boson Mass 

Since our models are defined in terms of weak-scale parameters, their implications for 
physical masses are typically clear, and we simply avoid regions of parameter space that are 
excluded by direct searches for new particles. 

The one exception to this rule is the lighter CP-even Higgs boson h. At tree level, 
m/j mz for tan/3 > 10 and rriA 3> mz-, and raising it above the LEP bound of 114.4 GeV 
requires significant one-loop corrections involving top squarks. In the minimal supersym- 
metric standard model (MSSM), for mi < 1 TeV, significant left-right top squark mixing is 
required to boost m^ to acceptable levels. In non-minimal models, however, it is possible 
to increase the lightest Higgs mass in other ways. For example, in the next-to-minimal 
supersymmetric standard model with a singlet scalar mass of Ms ~ 5/i, the Higgs boson 



mass may be as large as rrih ~ 125 GeV 74], even in cases where the MSSM contribution 



leaves m^ ~ mz- Although this introduces a new complex scalar and a neutralino that in 
general mixes with the rest of the neutralinos, the mass is sufficiently high that it effectively 
decouples. 

For these reasons, we set rrih = 114.4 GeV whenever the explicit one-loop corrections 
of our model are insufficient to drive it above that bound. Note that this redefinition is 
primarily important for sfermion masses disfavored by the EDM constraints discussed above. 
In addition, raising the Higgs boson mass generally decreases a , so simply removing models 
with too- low rrih would, if anything, raise the lower limit of our characteristic range of a . 



3. Extrapolation to High Scales 

Because these models are specified by weak-scale input parameters, there may be tensions 
associated with extrapolating them to high energies through renormalization group evolu- 



tion. Studies of non-unified Higgs models [60|, |63J have found that particular choices of fx 
and rriA produce an unstable Higgs potential at the high scale due to renormalization effects. 
This problem is generally avoided for fi ^ 1 TeV and large ttia- The exact requirement for 
rriA depends on other model parameters, but rriA ~ 1 TeV is generally large enough to avoid 
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this difficulty. A moderately large value of ttia ~ 300 — 500 GeV is also motivated by the de- 
sire for a sufficiently large relic density, with low values generically producing overly-efficient 
annihilation even without resonance effects. 

Large gluino masses, with M3 3> rriq, may also imply tachyonic squarks at the high 
scale. This is not necessarily a problem, as the low-energy theory could be located in a false 
vacuum that is stable on cosmological time scales, a possibility that has been considered in 
several studies (75-80[. Of course, all of these potential tensions assume a desert above the 
TeV scale and may be avoided by the introduction of new heavy fields without altering the 
phenomenological analysis presented here. 



III. RELIC DENSITY 

In this study, we require that the neutralino be a thermal relic, with Q x = 0.23. In this 
section, we analyze what this constraint implies for the gaugino and Higgsino content of 
the neutralino LSP, which will help us understand the results for neutralino scattering cross 
sections. 

The mass eigenstate x is a mixture of gaugino and Higgsino eigenstates given by 

X = a § B + a w W + a ku H u + a 6d H d . (13) 
With the assumption of gaugino mass unification, <C 1, since M 2 — 2Mi, and so x is 



primarily a Bino-Higgsino mixture. The exact forms for &jj £11*6 [81 



afj d _ \i (M 2 -m x ) (M 1 -m x ) - \m\ sin 2/3 \{M l -M 2 ) cos 2 6 W + M 2 - m x ] 



(14) 



ar. 



g mz (M 2 — m x ) sin 9w (/■* cos (3 + m x sin (3) 

3r„ -I- 7lf„ — m 1 

(15) 



a H u _ ~ m x (M 2 -m x ) (Mi-m x ) - m| cos 2 (3 [(M 1 -M 2 ) cos 2 6 W + M 2 - m x 



B 



mz (M 2 — m x ) sin 6\y (jj cos (3 + m x sin (3) 



For moderate to large tan/3, both sin 2/3 and cos 2 /3 are small, and the first terms in both 
numerators are the dominant contributions. This conclusion may be evaded only if /i <ti mz, 
which predicts a light Higgsino and is excluded by LEP bounds, or if Mi and m x are 
almost exactly degenerate, in which case x ~ B and \a^ \ 1. Thus, either x ~ B and 
\ag I <C 1, or % is a significant mixture of gaugino and Higgsino, and 



^w— (16) 



Given this fact, it is useful to define a single parameter, the "Higgsino-ness" 



which characterizes the size of both ag^ and a^ d when they are significant. Note that there 
is also an important dependence on the sign of /i, which determines the relative signs of Og 
and ajj ud . We will generally find solutions to Q x = 0.23 for both /i > and /i < 0, and so 
we will present results for both sign choices. 



S 



H 



W+ H 



H 



H+ 



H W- 
(a)Chargino-mediated 
Annihilation to W + W~ 

H 





h 



B W + 
(b)Higgs-mediated Annihilation 

to W+W- 





fR B 



h 



B Sl 

(c)Higgs-mediated Annihilation 
to Fermions 



* 
I 
I 
I 



h 



B h 

(d)Sfermion- mediated 
Annihilation to Fermions 



B JL 

(e)Sfermion- mediated 
Annihilation to Fermions with 
LR Mixing 



FIG. 1. Major processes contributing to S-wave XX annihilation. The diagrams are labeled with the 
appropriate gauge eigenstates of the neutralino to highlight the required B and H contributions. In 
|(d)| and (e) , additional diagrams exist with the exchange B f-> H and either _/x <-> or fi <-> f^. 



tempered neutralino scenarios |47|, |82 



Discussion of B-H mixing is familiar from focus point supersymmetry |40|, |42[ and well- 
However, we wish to stress that this mixing can be re- 



garded as the key parameter for all choices of parameters El, even in the other typically-defined 
"bulk" , "co-annihilation" , and "A-funnel" dark matter regions discussed in mSUGRA. Our 
approach, however, effectively separates this mixing from the correspondence of [i with scalar 
masses found in mSUGRA0, as well as avoiding other effects of running based on high-scale 
parameters which obfuscate the electroweak scale mass and mixing parameters. 

The dominant processes contributing to XX annihilation in most regions of parameter 
space are shown in Fig. [TlR The annihilation processes shown in Fig. 1(a) -Fig. 1(d) require 



a x with a H component and thus scale with ag. The remaining process of Fig. 1(e) does 



not require a H component, but it is generally suppressed relative to Fig. 1(d) because it 
requires left-right sfermion mixing. 



rn 



In Fig. [21 we show what Higgsino content is required to satisfy VL X 
t for fixed sign(/i), ttia, and tan f3. First consider Fig. 2(a) We see that 



0.23 as a function of 
when 



\ a H d \ 



they are significant, as discussed above. We also see the general trend that increases 
with increasing m x , because larger m x suppresses annihilation and thus ag must increase to 
compensate. This rise is interrupted three times: the dip centered at m x ~ 55 GeV is caused 

2 Even models without a unified gaugino mass, the W component of x is small so long as M2 is even 
moderately larger than Mi. 



GO 



G3] 



3 This is a general feature of NUHM scenarios 

4 Generally other i-channel processes, such as annihilation to ZZ, are sub-dominant. There are, however, 
specific regions of parameter space where resonant annihilation dominates for an intermediate particle of 
mass rn/ « 2m x or where co-annihilation is a major effect. We include all of these effects in our numerical 
results. 
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FIG. 2. Higgsino content of the dark matter particle x as a function of its mass. The Higgsino 
content is fixed by the requirement £l x = 0.23. Panel (a) shows ag, \ag |, and \ag | for fh = 2 TeV 
and both choices of sign(/u). The solid line corresponds to fj, > and the dotted line to fi < for 
ag, and the dashed lines correspond to /i > and dot-dashed lines to /U < for ag . Panel [(b)] 
shows for various values of fh for fi > 0. For both panels, we fix m/i = 4 TeV and tan /3 = 10. 



by the h resonance, where the annihilation cross section is enhanced kinematically, and so ag 
must decrease to compensate dynamically. In addition, there are two notable drops in ag at 
m x f=s m\y and m x pa m t , where annihilation is enhanced by the opening of a new annihilation 
channel. Finally, ag rises to a 1 at m x ~ 1.7 TeV. This corresponds to a maximum value of 
m x consistent with the relic density due to decreasing annihilation efficiency as m x increases, 
barring other effects (such as co-annihilation) raising the annihilation rate once more. There 
is also splitting between the cases of \i > and fi < 0, a result of the aforementioned 
dependence of ag on the sign of fi and the interplay of these parameters in annihilation. 

shows the dependence of on fh. For m x < mt the dependence on fh is 
weak, primarily because the process XX ~^ W + W~ is independent of fh and the sfermion- 
mediated processes are suppressed relative to this by Yukawa couplings. However, when 
XX ~^ is kinematically allowed, the t-mediated process may be similar in magnitude to 
the Higgs-mediated process. In this case, for fh = 2 TeV the Higgs-mediated annihilation 
is significantly larger than the squark contribution and causes a drop in across the ti 
threshold; for lighter fh the t-mediated contribution cancels the Higgs-mediated contribution 
to some degree, and thus less of a drop in ag is required. This effect is especially visible in the 
case of fh = 400 GeV, where there is very little effect from passing m x = m t corresponding 
to nearly exact cancellation between the two processes, and fh = 300 GeV, where ag again 
drops significantly because the t-mediated diagram is the dominant contribution. The curves 
for lower values of fh also show the effects of co-annihilation, with the value of ag dropping 
rapidly to a small value then remaining small up to the kinematic limit H The actual value 

5 These effects should also be present for the 1 TeV and 2 TeV curves, but were unreachable due to numerical 
instabilities in finding a valid window for the relic density at high m x . 



Figure 2(b) 
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may vary between zero and a finite but small maximum value based on the value of \i\ we 
display a value near the maximum to illustrate later effects upon the scattering cross section 
and to minimize the value of \x in the co- annihilation region at only ~ 2 TeV. 

It is instructive to explicitly relate the behavior of in Fig. [2] back to the fundamental 
parameters in the theory. Throughout most of the range of m x shown, 0.1 < < 0.9, 
which is associated with the focus point region where |Mx| ~ | /x | to allow for significant 
mixing. In particular, for that range of a^, 0.67 < \Mi/ fi\ < 1.05. The significant drops 
due to the Higgs resonance and co-annihilation correspond to a nearly pure B LSP with 
| Mi | 3> |/i|. We have set tua = 4 TeV, so there is no A-funnel region present in Fig. [2j 
but for a lower value of tua there would be a small region of significantly suppressed at 
m x ~ The presence of an A-funnel effectively obfuscates other interesting features 

in the behavior of and ultimately a , and this is part of our motivation for examining 
a large base value of throughout. The behavior for lower values of tua is discussed in 

Sec. ivra 



IV. SPIN-INDEPENDENT DIRECT DETECTION 
A. General Formalism 

The relevant interactions for spin-independent scattering are four-fermion neutralino- 



nucleon effective operators. Following closely the notation of Ref . [54| , these are 

C SI = \ N 4> x ip x 4> N i) N , (18) 

where N = p,n, and is determined by the underlying SUSY interactions that produce 
the operator. For a nucleus with charge Z and A total nucleons, this leads to a total cross 
section of 

a sl = ^ [X P Z + X n (A - Z)] 2 , (19) 

71 

where /i x = m x M^/ (m x + M^) is the reduced mass of the x-nucleon system. For neutralinos, 
typically \ p ~ A„, and the cross-sections add coherently, so the cross section can usually be 
reduced to 



a 



4/x 



2 



SI _ \2 A 2 



71 P 



XzA . (20) 



The dominant diagrams contributing to <r SI for mixed B — H are shown in Fig. [3j all 
of which have different coupling strengths for different quark flavors. Because of this, the 
cross-section has a strong dependence on the individual contributions of each quark flavor 
to the current ip q ip q in the nucleon. This quantity is parameterized as 

{N\m q ^ q \N) = f^M N , (21) 

where is the nucleon mass and is interpreted as the contribution of the quark q to 
the nucleon mass. 

Using this parameterization, the coupling of the neutralino to the nucleon becomes 

— = tf?-' (22) 
m N ^— j' y m q 
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(a)Higgs-mediated Scattering 




h 



Ql Qr 
(b)Squark- mediated Scattering 




X- 

ql qr 

Ql Qr 
(c)Squark-mediated Scattering 
with LR Mixing 



FIG. 3. Major processes contributing to a SI . The diagrams are labeled with the appropriate gauge 
eigenstates, and additional diagrams exist with the exchange B -H> H and either Jl *h> or 
/i o This exchange is independent in (a) and correlated in |(b)| and (c) 



where X q is the effective coupling of each individual quark to the neutralino through an 
operator analogous to Eq. f JT8|) above with iV — )■ q. For the light quarks u, d, and s, the 
values of can be taken from meson scattering experiments or lattice results. For the 
heavy quarks Q = c,b, t, the entire contribution derives from loop contributions through 
gluon exchange, and has the value [83 



v 




q=u,d,s 



(23) 



The parameters f^ d are reasonably well-known and small, but has been the subject of 
debate due to varying experimental and lattice results. In our analysis we take the values 
of f^ d from Ref. [54( and use the lattice result given in Ref. [84J for which gives 



f d = 0.033 



'd 

r-n 
Id 



0.042 



fp 

J u 
r-n 
J u 



0.023, /f 

: 0.018, /; 



0.05 
= 0.05 



(24) 
(25) 



These values imply /q = /q = 0.066, so the heavy quark contribution is quite important 
to the overall value of cr SI . We will discuss the impact of variations from these values in 
Sec. IVTEl 
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B. Neutralino-Quark Couplings 



The coefficients for spin-independent scattering of neutralinos off individual quarks, X q , 



have the form [85 

X u = 



2 (m|. - ml) 



m u a H u 

2mw sin p 1 



(it)* 



i2 



X 



2mw sin p 3 



# 2 m n 



4mvy sin (3 



(a-fir — tan 6 w a^) 



[a jj d cos a sin a + cos 2 a ) 



mf 



+ (—a^ d cos a sin a + ajj u sin 2 a ) 



n? 



H 



(26) 



2 (m 2 - 



(d)* 



2miy cos /9 



y »2 



X 



9 2 m d 
Amy/ cos /3 



(a^ — tan^iya^) 



mda ^d rr(<*) i j. /) rr (d) 

L/^ + - tan 9wa s U> 2 

2m w cos p 3 

(a fj u cos a sin a + a^ d sin 2 a ) — j 



+ (—CLjj u cos a sin a + a^ d cos 2 a ) 



m 



H 



(27) 



where a is the CP-even Higgs mixing angle, m qi and U^f are the squark masses and mixing 
matrix defined by 



M~ = U m 



g. „„„ „ v 

2 



m gi 
ml 



(28) 



and there is an implicit sum over i. The m^-dependent terms of Eqs. (1261) and (T5j 
arise from the diagram of Fig. 3(a), and the m^-dependent terms are from the diagrams 
of Fig. 3(b) and Fig. 3(c) The coupling for heavier quarks can be determined by simple 
mass replacement in the appropriate equation. These couplings can be altered somewhat by 
the inclusion of an additional CP- violating phase in the Higgs sector 86j, but we will not 
address this contribution in this work. 

We have assumed M| = m 2 l, which implies [7® = 1. It is instructive to consider the 
limit mjj ~ m^ — > oo. In this case, — 7r/2<0, and the couplings simplify to 



X u 



g 2 m u (ciyy — tan6*iya 



B) 



Amy/ sin (3 

g 2 m d (a w - tan 9 w a^) 
Amy/ cos (3 



an 



m 2 — m 2 

A. 



+ 



i jj sin j3 cos (3 + sin 2 (3 



inr. 



a H d 



m 2 — m 2 



+ 



i fy u sin (3 cos (3 + a & cos 2 (3 



mi 



(29) 
(30) 



The dependence of X q on the Higgsino-ness of the neutralino is manifest, with the coupling 
vanishing for a pure B or H and the maximum value generically achieved for a well-mixed 
case. Moreover, it is clear that the squark-mediated diagrams are sub-dominant relative to 
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the Higgs-mediated diagrams, except when ttv~— m x ~ m\, in which case their contributions 
become comparable. 

Assuming the Higgs-mediated contributions dominate, there is also clearly a dependence 
on the relative signs of ag and ag , both in the Higgs-mediated contribution alone and 
in the relative sign between the contributions. This contribution is either a suppression 
or enhancement in both types of quarks at once, with a characteristic suppression of the 
couplings when ag have the same sign relative to the values with opposite signs. Moreover, 
in the same-sign case there is the possibility of cancellation between the two contributions 
further reducing <r SI , while for the opposite-sign case there is an enhancement. Equation (Fl6|) 
implies that ag u and ag will have the same sign for \i < and opposite signs for \i > 0. 
This means there will be a characteristic enhancement in the \x > value of a sl over the 
H < value for a given set of other parameters. This holds true over most of our parameter 
space. The most significant exception is when the LSP is a nearly pure Bino and Eq. f fT6|) 
doesn't hold. In this case, < 1, which causes a severe suppression in \ q and thus in cr SI 
as well, regardless of the sign of /i. 

V. CHARACTERISTIC SPIN-INDEPENDENT CROSS SECTIONS 

We are now ready to present our main results for the spin-independent neutralino-nucleon 
cross sections in our model framework. In Fig. HI we present the predicted values for a as 
a function of m x for models with the correct thermal relic density Q x = 0.23. Results are 
given for various fh, both signs of /i, = 4 TeV, and tan (3 = 10. The dependence on 
and tan /3 is examined in Sec. IVI Al 

A. Characteristic Features 

Figure H] has a number of striking characteristic features. First, aside from the results for 
fh = 300 GeV and /i < 0, which we discuss below, the predicted cross sections a SI are almost 
independent of dark matter mass. For fh = 2 TeV, this may be understood as follows: in 
this case, sfermions are decoupled, and scattering occurs through the light Higgs boson. 
These cross sections scale as ~ a 2 g/m\, but the Higgsino-ness is typically constrained by the 
thermal relic density to the range 0.3 < ag < 0.8, as shown in Fig. El and the Higgs boson 
mass is constrained to the narrow range 114.4 GeV < m/j < 130 GeV, given the current 
experimental bound. As a result, for m = 2 TeV, = 4 TeV and tan/3 = 10, <r SI is in the 
narrow range 2 zb < <j si < 20 zb for all m x > 70 GeV. 

Perhaps more surprising is that this behavior persists to much lower sfermion masses. 
As evident in Fig. HJ even for m = 400 GeV, the characteristic range of <r SI is extended 
only to 1 zb < <t si < 40 zb for m x > 70 GeV. Such low values of m are not far from 
current bounds on squark masses from direct searches at the Tevatron [87]. In this sense, 
current bounds from colliders already imply that sfermions are essentially decoupled from 
neutralino phenomenology, except in fine-tuned cases discussed below, greatly simplifying 
SUSY parameter space and focusing predictions. 

In more detail, we also see that the a sl curves of Fig. H] show a generic similarity to the 
afj curves of Fig. [2] as suggested by the form of the coefficients X q . There is a peak and 
subsequent drop-off in the distribution at m x ~ mw and again at m x ~ m t corresponding 
to the fluctuations in the ag, and most of the curves demonstrate some rise after m x ~ m t . 
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FIG. 4. Spin-independent neutralino-nucleon cross sections. The cross sections a for models with 
thermal relic density f2 x = 0.23 as functions of neutralino mass m x for the various sfermion masses 
indicated, \x > (solid) and fi < (dashed), tua = 4 TeV and tan/3 = 10. The shaded region 
between the solid and dashed curves are obtained for complex values of fi that interpolate between 
the two real values (see Sec. IVID|) . 

For rh = 2 TeV, a SI drops of at m x < 1.7 TeV as the neutralino becomes nearly pure H. 
This because, in the absence of left-right mixing, the dominant scattering processes require 
X to have both significant B and H components as shown in Fig. [3j 

Finally, for the reasons discussed in Sec. IIVBI the value of a for /i < is almost always 
suppressed relative to that for fi > 0, usually by at least a factor of two. This effect is 
enhanced for lower m. This follows directly from the form of X q above, with decreasing 
scalar masses bringing greater enhancement for fi > and greater suppression for \x < 0. 
Although this relation is not strictly true, due to the more pronounced dip in cr SI | M<0 for 
larger m at m x ~ m t , the suppression in a for rh = 2 TeV is relatively minor. 

B. Fine-tuned Cancellations 

The cross section er SI may move outside its characteristic range in certain fine-tuned 
regions. Given our simple, phenomenological model framework, we may identify the fine- 
tuned possibilities relatively easily. 

First, for 50 GeV < m x < 65 GeV, <r SI may be suppressed by orders of magnitude, 
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consistent with the suppression of shown in Fig. |2j This behavior is visible at the left- 
hand edge of Fig. HI For these masses, rrih rs 2m x to within a few GeV, XX annihilation is 
enhanced by the Higgs boson resonance, and so the correct thermal relic density is obtained 
even for Bino-like neutralinos. In this fine-tuned case, <r SI may be so small that neutralinos 
will escape all direct detection searches for the foreseeable future. Note, however, that this 
possibility is disfavored not only by fine-tuning, but will soon be probed by chargino and 
gluino searches, since these are sensitive to LSP masses around 50 GeV under the assumption 
of gaugino mass unification. 

The regions in Fig.HJmust, of course, be cut off at m x = fh; for m x > fh, the neutralino is 
not the LSP. In fact, the behavior of a sl changes significantly for values ~ 20 — 50 GeV below 
fh. The reason for this is that the annihilation is enhanced for m x ~ fh by co-annihilation 
effects, quickly suppressing the relic density once co- annihilation becomes important. Al- 



though aft may compensate for this effect, the ag curves in Fig. 2(b) demonstrate rapid 
change as m x approaches the end of the curve, indicating that x wm quickly become a 
pure B or H to suppress annihilation somewhat before m x = rh. After this point the cor- 
rect relic density is unattainable. However, in contrast to typical assumptions regarding 
co-annihilation, the value of <r SI may actually increases greatly approaching the kinematic 
limit. This is due to the unified sfermion mass, which results in co- annihilation with squarks 
as well as sleptons, and thus there is a near- resonance enhancement of <r SI if remains 
small but constant. This effect is largely illustrative, as most realistic scenarios predict slep- 
ton masses at least somewhat smaller than squark masses, and the exact values of cr SI are 
undetermined as the exact value of ag is not set in the co-annihilation region. 

For rh = 1 TeV, 500 GeV, and 400 GeV, there is a small region with a significantly 
suppressed value of cx SI in the \i < curves for m x just below the co-annihilation cutoff. 
To understand this feature, consider the form of A^ from Eq. (T291 . For large tan/3, sin a is 
small and negative, so A^ = for 



2 a H d 2 o -in 2 , . > , , 

m 7 — m ~ - m h ~ m h ~ 10m h (31) 

x ag cos a sin a m x sm2a 



for characteristic values of a and fi/m x . For even moderately larger values of rh, the re- 
mediated contribution is sub-dominant until m x becomes sufficiently large. For m 2 x ~ m 2 - — 
10m|, X d vanishes at a specific value then flips sign for larger values of m x . For X u = 0, the 
corresponding condition is 

1 

m l ~ m l ~ o a : m l ; ( 32 ) 

A cos z a — — cos a sin a 

m x 

so X u cannot vanish for m x < fh. The highly suppressed region is the result of cancellation 
between contributions from up-type and down-type quarks in this region, the result of the 
squark-mediated contributions to <r SI in Fig. |3]becoming similar in size to the Higgs-mediated 
contribution for m x ~ fh. 

In fact, this effect explains the suppression of <r SI for fh = 300 GeV across a large range 
of m x . In this case fh 2 ~ 10m^ regardless of the value of m x , so the Higgs- and squark- 
mediated contributions are always comparable. This leads to a general suppression in <j SI 
until m x ~ fh, where a sl actually increases to the characteristic value because the squark- 
mediated contribution dominates. 
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FIG. 5. Theoretical predictions for o~ SI and experimental bounds. The shaded regions are as in 
Fig. [U along with various current and project bounds. The current upper limit from CDMS [f| 
skirts our projected cross sections, and recent XENON100 10[] results have begun to constrain our 
characteristic region. Most or all of the characteristic regions will be probed by the projected near 
future sensitivities of XENON100 [ll|, XMASS 0], LUX (s^S, and superCDMS with 100 kg 
fiducial mass at SNOLAB [13]. 



C. Prospects for Discovery 



The primary implication of Fig. H] is that simple SUSY models typically predict 1 zb < 
cr SI ^ 40 zb. This prediction becomes even more narrow for large sfermion masses, which are 
required for 0(1) phases to be consistent with the EDM constraints discussed in Sec. Ill B 1[ 



The reason lies in the decoupling of squark- mediated processes in Fig. 3(b) and Fig. 3(c 



leaving only the Higgs mediated process in Fig. 3 (a 



Such cross sections are just below current bounds for most masses, and recent XENON100 [10 
results have begun probing our characteristic region. These predictions along with current 
experimental bounds and expected near future sensitivities, are shown in Fig. [5j We see 
that the expected range is just below current limits, but the characteristic range will be 
probed by experiments that reach the zeptobarn scale in the near future. 
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FIG. 6. Dependence of a SI on Higgs potential parameters. The cross section o~ sl is given for various 
niA and tan /3 = 10 in (a) and various tan/3 and toa = 4 TeV in |(b)| In all cases, m = 2 TeV, and 



the upper contour is for [i > and the lower contour is for \x < 0. 



VI. DEPENDENCE ON MODEL PARAMETERS AND ASSUMPTIONS 

In this section we examine the effects of relaxing several of our particle physics model 
assumptions, varying the strange quark content of the nucleon, and modifying astrophysical 
inputs that enter our analysis. 



A. Variations in tua and tan/3 

As discussed in Sec. |ITJ our Higgs potential is parameterized by /i, m^, and tan/3. After 
imposing the constraint Q x = 0.23, the magnitude of /i is determined, and we have discussed 
the dependence of our results on the choice of sign(/i). Figure [6] shows the effects of varying 
the remaining input parameters of the Higgs potential: ttia and tan /3. 

As noted in Sec. Ill B 3\ overly efficient annihilation and stability of the Higgs potential at 
high scales generally motivate at least partially decoupled values of niA ~ 300 — 500 GeV. 
Up to this point in our analysis, we have considered the complete decoupling limit with 
rriA = 4 TeV. Here we consider the effects of a lower value of tua by showing results for 



tua = 1 TeV in Fig. 6(a) The primary effect is a strong suppression of <r SI for m x ~ 
For such values of m x , XX annihilation is enhanced by the S"-wave A resonance. This effect 
is similar to the light Higgs resonance discussed previously, as it forces x t° De nearly pure 
B to compensate for the strong annihilation enhancement and thus suppressed a . This 
behavior is present in the A-funnel region at large tan (3 in mSUGRA, but we see here that it 
is generically possible for all values of tan /3 and the position in parameter space is generically 
unconnected to the sfermion masses. There is also a slightly less pronounced resonance 
involving the heavier CP-even Higgs which produces the less distinct second downward 
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spike for the /i < curve. The value of <j si drops somewhat for m x > 1 TeV due to the 
availablitiy of heavy Higgs final states in annihilation. The availability of AA final states 
also increases the maximum value of m x by increasing annihilation efficiency Finally, there 
is a slight broadening of the values of <r SI due to the increased contribution from scattering 
through the heavy CP-even Higgs; this effect will become more pronounced for even lower 
rriA, but the aforementioned lower limit on prevents this effect from causing serious 
suppression. 

For variation in tan /3, the characteristic range of values for cr SI narrows for larger values 
of tan j3. This can be understood as follows: for the neutralino mass matrix, in the limit 
tan/3 — > oo, the sign of [i may be switched by the field redefinition Hd — > —Hd without 
affecting other matrix entries. We therefore expect the fi > and fi < predictions to 



become similar for large tan /3. Conversely, as shown in Fig. 6(b) , the range of cr SI becomes 
wider for lower values of tan /3; for tan /3 = 5, roughly the minimum value allowed by null 
results from Higgs boson searches, the lowest value of a is reduced by a factor ~ 2 from 
the tan (3 = 10 results. 

For very large tan /3 and low value of ttia, there is a possibility of a general suppression 
due to cancellation between light and heavy Higgs contributions to cr SI . For tan/3 — » oo, 
sin a — > and cos a — > 1. Then, if the squark mediated contributions are neglected, the 
couplings in Eqs. ( 126]) and ( p7|) become 

2 

A " * ~^lkp ^ ~ tan6wa ^ (33) 

Xd ^^^^' t ^ 6wa ^ H - (34) 

Then, under the (very) rough assumption of equal portion of up- and down-type quarks in 
the nucleon, there is a cancellation between these contributions when X u /m u + Xd/rrid — 0, 
which gives 



!^~ 4 /°5i tan/3 tan/3 (35) 

This allows for a possible cancellation for [i < with little dependence on m x . For tan/3 = 10 
the required value is m# ~ ~ 300 — 400 GeV, which strains the relic density calculation; 
however, for tan /3 = 50, the cancellation occurs at larger values of ~ 700 — 800 GeV. 
However, even for large tan /3 significant suppression only occurs at very specific values of 
m,A, as for lower values of rriA the contribution to scattering from the heavy Higgs dominates 
the cross section. 



B. Un-unified Sfermion Masses 

Although convenient for studying the general behavior of <r SI , our assumption of a unified 
weak-scale sfermion mass is not motivated by a connection to a high-energy theory. Indeed, 
most high energy theories predict < rrig, either resulting from a unified scalar mass at 
the GUT scale, as in mSUGRA, or as a remnant of sfermion masses generated by gauge 
dynamics, as in gauge-mediated supersymmetry breaking models or anomaly-mediated su- 
persymmetry breaking models with a unified scalar mass. 



To explore the consequences of lighter sleptons, in Fig. 7(a) , we compare our standard 
scenario with rh — 1 TeV to a scenario with a unified weak-scale squark mass nig = 1 TeV 
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FIG. 7. Dependence of o~ SI on the unified sfermion mass assumption. Cross sections cr SI are given 
for fj, > (upper) and fi < (lower), tua = 4 TeV, and tan/3 = 10. Sfermions masses that are not 
specified are set to rh = 1 TeV. 



and a separate unified weak-scale slepton mass = 200 GeV. The lower value for 
reduces the value of a^ 1 across the entire range of masses by ~ 20 — 30%. This is due to 
the effect of t-channel ^-mediated diagrams in \ annihilation reducing ag and providing no 
additional to cr SI . As expected, the spike of high <r SI associated with squark co- annihilation 
in Fig. H] figures is absent. While it is possible that even smaller values of m~ t could produce 
much stronger suppression of <r SI over a large range of m x (that is, not just in co-annihilation 
regions), the value of = 200 GeV is already in significant tension with electron EDM 
constraints. Thus we conclude that the main effect of lower is to reduce the range of 
m x for which the \ can be the LSP, except perhaps for very low values m~ t which imply 
significant fine-tuning in CP phases. 

Another well-motivated non-unified sfermion mass possibility is to have heavy first and 
second generation sfermions and light third generation sfermions. Precision measurements 
strongly constrain the mixing in the first two generations, but the mixing between the first 
two generations and the third generation is less well-constrained. An interesting possibility, 
then, is that the squarks of the first two generations are heavy to limit contributions to 
low-energy observables through decoupling, while the third generation squark masses are 
light to avoid fine-tuning in the electroweak scale [90l - The impact of such a sfermion 



mass pattern is shown in Fig. 7(b) , which contrasts results for rh = 1 TeV and the case 
where m^i = 500 GeV, but all other sfermions remain at 1 TeV. Beyond lowering the upper 
limit on m x such that x remains the LSP, a lower value for mi reduces the suppression of 
cr SI associated with ti annihilation, as discussed in Sec. [TV] There is also a slight widening 
of allowed values of <r SI , as do contribute to a sl , but this widening is not as significant 
as that caused by lowering rh. Overall, the effect is larger than the effect of varying slepton 
masses, but again, it is not very significant. 
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FIG. 8. Dependence of a on left-right sfermion mixing. Cross sections a are given for no left- 
right mixing and maximal left-right mixing for 3rd generation sfermions, fi > (upper) and fi < 
(lower), rh = 2 TeV, uia = 4 TeV, and tan (3 = 10. 



C. Left-Right Sfermion Mixing 



Although we have assumed negligible left-right sfermion mixing, there may in general be 



some mixing, and this will alter both \X annihilation through the process of Fig. 1(e) and 



(j SI through the process of Fig. 3(c) Even if left-right mixing is only significant for third- 
generation sfermions, as is typically true, in principle this may affect <t si by modifying XX 
annihilation in the early Universe or by changing the heavy quark contributions to scattering 
off nucleons. 

Figure |8j shows the effect on cr SI of left-right mixing in third generation sfermions. The 
mixing here is parameterized by the angle 9q, and the two cases shown are those with no 
mixing (sin 29 q = 0) and maximal mixing (sin 29q = 1), for which each mass eigenstate has 
equal parts /t and These are achieved by fixing the A t;b)T parameters appropriately. 
Although the latter case does produce some splitting between the / mass eigenstates, for 
the degenerate masses considered here, a highly mixed scenario can be achieved without a 
kinematically significant splitting for our analysis!! 

As seen in Fig. [BJ for rh = 2 TeV, left-right mixing may enhance cr SI by ~ 15%. The 
effects of left-right mixing on the relic density calculation are minimal, but contributions 
from Bq^^R — > BqR t L and HqL,n —> Hqn,L produce an enhancement in cr SI directly. The 
effect may be larger for lower rh, but in any case, we find that the overall effect does not 
significantly impact conclusions about the characteristic range of a SI . The value of cx SI is 



enhanced as \ becomes pure H at high mass by a process similar to Fig. 3(c) with H initial 
and final states, but this process is yukawa suppressed and <r SI still drops below 1 zb in the 
highly H region. 

6 The masses are not exactly degenerate, given corrections from D-terms, but even including these, a pair 
of completely mixed states can be created with a mass splitting of < 2 — 3 GeV for the values of m 
considered. These corrections are kept even in the no mixing case to prevent unintentional left-right 
mixing from numerical effects. 
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D. Explicit Inclusion of CP- Violation 



Although the desire for 0(1) CP-violating phases motivates large sfermion masses, we 
have performed our analysis is a manifestly CP-conserving framework. This is a cause for 
some concern, as CP violation could impact relic density and scattering calculations, which 
could, in principle, suppress <r SI . In our model framework all CP violation is found in the 
phases of the parameters Mi and fi. We adopt a convention where all CP violation is 
contained in the phase of the \x parameter 

The effects of general phases on both Q and cr SI were discussed by the authors of Ref. j§5|, 
who found that phases can significantly enhance the relic density and suppress <r SI . (See 
also Ref. (96[.) For a highly- mixed neutralino LSP like the case considered here, varying 
the phase causes a suppression in both annihilation and scattering cross sections, leading 
to an enhancement in the relic density by a factor of ~ 5 and a suppression in cx SI by a 
factor of ~ 50. However, in that analysis the fl x and <r SI were treated as two independent 
parameters. In our analysis, with an appropriate relic density applied as a constraint, these 
effects tend to cancel each other: to achieve a viable relic density, ajj must be increased to 
enhance annihilation, which will in turn enhance <r SI . 

To determine the effect, we note that the relic density is nearly minimal for 0^ = and 
maximal for = ir, with intermediate values for (p^ smoothly interpolating between these 
two limits. The value of a has opposite behavior, with a maximal value for 0^ = and 
minimal value for = ir. The two limiting cases are, then, precisely the cases \x > and 
< that we have considered throughout our analysis. Thus the inclusion of an arbitrary 
CP-violating parameter should result in a simple interpolation between the cases of fi > 
and n < as alluded to in Sec. [V] without resulting in values of <r SI outside of this range, 
except perhaps in finely-tuned portions of parameter space. 

E. Strange Quark Content of the Nucleon 

As discussed in Sec. HV] the value of cr SI is strongly dependent on the individual quark 
contributions to the nucleon, f^. Of these, the quantity that has traditionally had the 
most uncertainty is f^, with quoted values ranging from 0.05 to 0.23. Recent lattice results 
suggest the lower value, and we have used this value throughout, but it is useful to examine 
the variation of cr SI as is varied. Throughout this section we use the assumption /| ps 

To understand the general dependence on f^, consider the quantity 



which can be interpreted as the total quark contribution to the nucleon mass. The (perhaps 
unreasonably) high value = 0.23 implies f^ t ~ 0.44, and the low value = 0.05 implies 
/tat — 0.30. Since fg is largely unchanged as decreases further, even the unrealistic 
case of ff = would still have / t ^ t 0.25. Thus, without considering further possible 
cancellations, the difference in <r SI between between no strange quark contribution to the 
nucleon and the largest proposed value should be about (0.44/0.25) 2 ~ 3. A more careful 
variation is shown in Fig. |9] for m = 2 TeV and 300 GeV. The 2 TeV case matches the 
sketched analysis above quite well — the value of cr SI for = 0.05 is suppressed by ~ 2 
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FIG. 9. Dependence of a SI on the strange quark content of the nucleon. Cross sections are given 
for the different values of /f indicated, = 4 TeV, and tan j3 = 10. In panel (a) , rh = 2 TeV 
and \x > (upper) and /i < (lower). In panel [(b)| m = 300 GeV, and only the fj, < curves are 



shown, as the ^ > curves have much the same behavior as in panel (a) 



relative to the value for / s = 0.2, with the value for =0.1 lying in between, in good 
agreement with the prediction based on f^ t . Note that, for rh = 2 TeV, our default value 
of yields the lowest a , and other possible values imply improved prospects for direct 
detection. 

For rh = 300 GeV the suppression for low f^ is much more pronounced, at least in the 
case of n < 0. As discussed in Sec. HVl for // < the sfermion contribution to scattering is 
similar in size to the Higgs contribution, opening the possibility of cancellation regardless 
of the value of m x . In fact, this also opens up the possibility for additional cancellations 
between X q for different quarks based on relative values of f%. The cancellation is most 
prominent when f^ is comparable to the contributions of other quarks; for larger values 
of the strange quark is the dominant contribution to nucleon mass, which raises the 
generic value of <r SI . Even for larger f^, however, <r SI for rh = 300 GeV remains suppressed 
relative to higher values of rh due to the squark-Higgs contribution cancellations discussed 
in Sec. El 



F. Variations in Relic Density 

It is possible that neutralinos form some, but not all, of the dark matter. One may 
consider models with Q x < 0.23, which may have larger o . Of course the direct detection 
signal is proportional not to a SI , but to Q x a SI . In fact, as is well-known, for thermal relics, 
the relic density is inversely proportional to (cr^v) , the thermally-averaged annihilation cross 
section. If a SI and <T an scale together, then the gain in a SI is exactly compensated by the 
reduction in fl x . 
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FIG. 10. Dependence of a on Q x . Cross sections a are given for the values of Q x indicated, 
fjL > (upper) and fi < (lower), rh = 2 TeV, niA = 4 TeV, and tan/3 = 10. 



In Fig. HOI we explore the impact of variations in Q x . For concreteness, we consider 



models with Q x = 0.184. As expected, in Fig. 10(a), we see that these models have larger 



rSI 



a"'. However, in Fig. 10(b) , we see that this is almost exactly compensated by the reduction 
in Q x except at high m x , leaving the direct detection signal invariant. At very high m x , 
the neutralino becomes nearly pure H and reaches its maximum annihilation efficiency 
more quickly for the lower value of Q. For relatively small deviations from the assumption 
Q x = 0.23, then, the direct detection signal is not impacted by the exact value of fl x assumed 
except for large m x . This result may not hold for large deviations when x constitutes a very 
small portion of the dark matter, but in this case the theoretical motivation for \ dark 
matter is somewhat strained. 



G. Galactic Small-Scale Structure 



The exclusion limits on cr SI are also affected by the local density of DM. Although a 
naive approximation suggests a spherical DM halo with uniform density at a given radius, 
small scale structure exists, and leads to a few regions with higher than average density, 
compensated for by a lower than average density in the rest of the space. This has been 



studied in Ref. [97] . which concluded that while the reduction in general number density 
from clumping might be as much as an order of magnitude, the more likely possibility is a 
reduction by less than a factor of two. This implies that experimental limits are subject to 
astrophysical uncertainties, but does not significantly alter our conclusion that near future 
experiments will probe the heart of neutralino parameter space. 
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H. Extension to Non-minimal Scenarios 



Because we work in a manifestly low-energy framework, our results can easily be gen- 
eralized beyond the MSSM to models that include new superfields with TeV-scale masses. 
Although such fields will in general affects the annihilation and scattering of \ dark matter, 
we argue that a sl remains generically in the range 1 zb ^5 0" SI < 40 zb. 

To illustrate this point, consider first the impact of a new TeV-scale particle X on anni- 
hilation. Assume that X couples to x an d some set of other MSSM particles with a coupling 
of order aw The effect of X will be most pronounced if there is resonant annihilation 
when mx ~ 2m x or if there is significant co-annihilation when m>x ~ rn x . Both cases force 
a severe suppression in to compensate and produce a sufficiently large relic density^, 
in turn causing a large suppression in <r SI similar to that seen in the h and A resonances 
discussed above. However, outside of these fine-tuned regions of parameter space, new X- 
mediated processes must be roughly the same size as existing processes to be relevant and 
will generally decrease (or possibly increase in the case of cancellation) ag relative to the 
MSSM case. Although this will produce a suppression (or possibly enhancement) of cx SI , the 
additional effect should not generically be greater than the factor of 2 — 3 that is caused by 
the introduction of W + W~ or tt annihilation. 

In scattering, if X couples to quarks, it may induce additional contributions to a SI . 
However, although the mass and interactions of X may be chosen to significantly suppress 
cr SI , such an effect requires a fine-tuning of the mass and couplings of X to produce a 
cancellation similar to the case of m = 300 GeV discussed in Sec. [TV] Indeed, the low value of 
m x required to cause a severe suppression directly through scattering effects should naively 
interfere with low-energy observables due to its couplings to quarks, and such operators 
are constrained by collider searches [29-3l||. Even if constraints on the mass of X can be 



avoided by tuning couplings, a sufficiently light mass for X places it generically within 
detection reach of the LHC. 



VII. SPIN-DEPENDENT CROSS SECTIONS 

We have focused in this work on spin-independent direct detection, as this provides 
bright prospects for the discovery of neutralino dark matter. It is interesting, nonetheless, 
to examine the prospects for spin-dependent cross sections, which is probed both by direct 
detection experiments, and by indirect detection searches for neutrinos. 

Spin-dependent scattering proceeds through the effective operator 

ociXl^XQa^Qi ■ (37) 

The structure of these operators implies that the contributing diagrams will be similar to 
those shown in Fig. [3J except with Z-mediated scattering replacing Higgs-mediated scatter- 
ing, and the replacement B -H- H and associated fermion L <-)■ R replacements on one leg in 



7 If reducing does not compensate for the resonant annihilation, then that region will generally not 
provide a viable relic density. 
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both Fig. 3(b) and Fig. 3(c) The associated coefficients have the form 85 
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(39) 



with the same notation used in Sec. IIVBI 

The spin- dependent couplings will exhibit a similar suppression of sfermion-mediated 
diagrams relative to Z-mediated diagrams as seen in the spin-independent case, with the 
smaller value of mz relative to m^ leading to the naive expectation that the squark mass 
necessary to cause serious suppression should be lower than the 300 GeV mass necessary 
to suppress cr SI . Such a small mass is ruled out by collider bounds. Moreover, as long as 
afj are non- negligible, the squark-mediated contribution to the cross section suffers an 
additional suppression relative to the Z-mediated contribution of 



sin 4 # 



— \a 



(40) 



in contrast to the minor suppression of the Higgs-mediated contribution relative to the 
squark-mediated contribution due to the mixing angle a for a SI . 

Predictions for <r SD for neutralino-proton scattering are given in Fig. [HI using the same 
parameters as in Sec. |V] and three values of m. In all three cases, the dominant effect is 
a drop in cr SD with increasing m x . This is due to the proportionality of the Z-mediated 
contribution to the factor 



— la 



as demonstrated in Fig. 2(a), the values of \a 



converge with increasing m x , thus suppressing <t sd . The dependence on m is mild until 
higher values of m x where the Z-mediated contribution is self-suppressed to the level of the 
sfermion-mediated contribution, and even the differences in <r SD at m x < 100 GeV are due 



to the slight variation in with mass shown in Fig. 2(b) rather than cancellations between 
scattering processes. 



The limits plotted in Fig. [TT] correspond to upper limits on <t sd from the AMANDA 98 
and IceCube [99] experiments, along with projected upper limits from a 10 year run of 
IceCube with DeepCore, an extension of IceCube with tighter string and optical module 
spacing jlOCrf ] and a 4 kg COUPP, assumed to run background free for 3 months at a deep 
underground site jlOl ] . The entire characteristic region plotted lies outside current upper 
limits, but the projected IceCube limits intersect the allowed region for 60 GeV ^5 rn x ~ 
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FIG. 11. Characteristic values for spin- dependent neutralino-proton scattering cross sections o~ SD . 
The shaded regions are cr SD for the values fh indicated, with /x > (solid border) and /i < (dashed 



border), txia = 4 TeV, and tan/3 = 10. Also shown are current limits from the AMANDA 98] 
and IceCube [9^] experiments, and the projected sensitivities of a 10- year run of IceCube with 



DeepCore 100[] and a 4kg COUPP running at a deep underground site 10l| 



200 GeV for most values of fh, and the fh = 500 GeV region may be within IceCube's 
sensitivity region up to m x ~ 300 GeV. 

We conclude that, at least for the neutralino DM models considered here, spin-dependent 
searches are indeed not as promising as spin-independent searches in the near future. How- 
ever, for some parts of parameter space, spin-dependent searches may also see a signal in 
the coming years. Of course, if a signal is seen in spin-independent detection, confirmation 
through spin-dependent signals will be an essential supplement. 



VIII. CONCLUSIONS 



Experimental sensitivities to the scattering of dark matter off normal matter are improv- 
ing rapidly. In this study, we have addressed the question of what this progress implies 
for the discovery prospects for neutralino dark matter. To determine characteristic spin- 
independent and spin-dependent neutralino-nucleon cross sections, we have worked not with 
specific high-energy models, but rather in a phenomenological framework incorporating only 
the basic motivations for SUSY (the gauge hierarchy problem, force unification, and the 
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WIMP miracle) and experimental constraints (direct searches and bounds on flavor and CP 
violation) . 

SUSY parameter space is large and complicated, but our strategy has been to start in a 
simple corner and work out from there. The conclusions of this study may be summarized 
as follows. We assume gaugino mass unification, as motivated by gauge coupling unification, 
and, to begin with, also heavy squarks and sleptons, as motivated by stringent constraints 
on flavor and CP violation. All of neutralino dark matter phenomenology then depends 
only on the Bino and Higgsino mass parameters Mi and /x, and the heavy Higgs scalar mass 
rriA and tan (3. Away from Higgs resonances, the dependence on the latter two parameters 
is weak. Requiring the correct thermal relic density removes another degree of freedom, 
leaving only one, which may be taken to be the neutralino mass m x . For m x > 70 GeV, the 
spin-independent cross section lies in the range 1 zb ^ a SI ^ 40 zb and is highly insensitive 
to m x for the reasons detailed above. The assumption of a unified sfermion mass disregards 
the complicated spectrum details associated with any high-scale model, but according to 
the discussion above such effects should be sub-dominant in any case except in very special 
cases such as co-annihilation. While the mass splittings of full spectra would certainly alter 
a specific prediction for a SI , such effects will not alter the characteristic value of <r SI studied 
here. 

From this result, we then generalize to other regions of parameter space. Surprisingly, the 
upper and lower limits of the characteristic range only vary by factors of ~ 2 for sfermion 
masses as low as m ~ 400 GeV. For very light sfermions near their experimental bound 
m ~ 300 GeV, large cancellations between contributions from up- and down-type quarks are 
possible. We have also identified and analyzed two other fine-tuned scenarios with greatly 
reduced a sl , where annihilation is enhanced either by resonances when 2m x ~ or ttia or 
by co-annihilation when m x ph m. In these special cases, cr SI falls outside our characteristic 
range. 

We next considered the effects of un-unifying the slepton and squark masses, un-unifying 
the third generation squark masses, including significant left-right sfermion mixing, and 
including CP-violating phases in the relic density and <r SI calculations. In each of these cases, 
for fixed Q x , we find typically small variations in the allowed range of cr SI . We have also 
explored variations of the strange quark content of the nucleon and found that our assumed 
value of fg = 0.05 is conservative, with other values enhancing <r SI for m <^ 400 GeV. 
Finally, we also considered the possibility of mult i- component dark matter where neutralinos 
are just one significant component. In this case, neutralino annihilation must be enhanced to 
reduce the thermal relic density, but this also enhances a . These effects effectively cancel, 
leaving the experimental direct detection signal, which is proportional to Q x a sl , invariant 
and in our characteristic range. 

These results for cr SI are compared to experiment in Fig. [5J If dark matter is composed of 
thermal relic neutralinos, the prospects for near future experiments that are sensitive to the 
zeptobarn scale are truly promising. We have also analyzed the prospects for spin- dependent 
scattering. These are summarized in Fig. [HJ Although not as promising as spin-independent 
scattering, observable signals exist for dark matter masses m x < 500 GeV. 

Our results are complementary to those based on high-energy frameworks. In the con- 
text of mSUGRA, for example, the fine-tuned possibilities of resonant annihilation and 
co-annihilation also exist and are known to predict low a sl . This results of this study show 
that these phenomena are decoupled from other SUSY phenomenology typically related to 
them in mSUGRA. For example, the A resonance may be realized at any tan/3, as opposed 
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to only high tan/3 in mSUGRA. Nevertheless, these possibilities remain fine-tuned, in the 
sense of requiring specific relations between superpartner masses that are not suggested by 
experimental constraints or foundational motivations of SUSY. 

At the same time, this study shows that results from the focus point region of mSUGRA 
are robust. As is well-known, in the focus point region, the neutralino is a well- mixed Bino- 
Higgsino combination and, for this reason, rates for both direct and indirect detection are 
large, with predictions for spin-independent cross sections in the range a sl ~ 1 — 10 zb. The 
results of this study show that similar cross sections are achieved in a more general setting 
throughout parameter space, largely independent of assumptions about GUT boundary con- 
ditions, A parameters, tan/3, or even the details of the slepton and squark spectra, provided 
fh ^ 400 GeV. Large direct detection signals and light sfermions are therefore perfectly con- 
sistent, with positive implications for both dark matter experiments and colliders. Given the 
generality of our model framework, the results for characteristic <r SI will apply to a diverse 
array of high-energy models. 

In conclusion, the field is reaching a critical juncture in the search for neutralino dark 
matter. We have identified a range for characteristic cross sections for neutralino-nucleon 
scattering. These are by no means a lower limit on possible value of cr SI ; indeed several 
possibilities for much smaller cross-sections were identified in this work. Nevertheless, even 
in the most pessimistic scenario, the identification of a characteristic cross section range has 
important consequences. If no signal is seen at the zeptobarn scale, then a model with ther- 
mal relic neutralinos must contain some amount of fine-tuning. This could be accomplished 
by enhanced annihilation through resonances or co-annihilation, as noted above. Alterna- 
tively, one could introduce new particles that induce neutralino-quark couplings that cancel 
existing contributions to cr SI . As we have seen, however, the scale of existing contributions is 
set by the light Higgs boson mass rrih, not by a superpartner mass. These new contributions 
would typically be associated with light, colored fields, then, as in the case of the ~ 300 GeV 
squarks we discussed in detail, and would be within reach of collider experiments. 

Of course, the most straightforward and natural conclusion is that if thermal relic neu- 
tralinos are a significant component of dark matter, their spin-independent scattering cross 
section likely lies in the range 1 zb < <x SI < 40 zb, just beyond current sensitivities. Current 
and near future experiments will therefore see a signal as they probe down to the zeptobarn 
scale, and the next generation of experiments will compile the statistics required to defini- 
tively establish the signal on different targets, constrain the dark matter's mass, see signals 
in directional detectors, and usher in the era of dark matter astronomy. 
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